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Abstract 

A new understanding of the consequences of how ecosystem elements are interconnected is emerging from 
the development and application of Ecological Network Analysis. The relative importance of indirect effects 
is central to this understanding, and the ratio of indirect flow to direct flow (I/D) is one indicator of their 
importance. Two methods have been proposed for calculating this indicator. The unit approach shows what 
would happen if each system member had a unit input or output, while the realized technique determines 
the ratio using the observed system inputs or outputs. When using the unit method, the input oriented and 
output oriented ratios can be different, potentially leading to conflicting results. However, we show that the 
input and output oriented I/D ratios are identical using the realized method when the system is at steady 
state. This work is a step in the maturation of Ecological Network Analysis that will let it be more readily 
testable empirically and ultimately more useful for environmental assessment and management. 

Keywords: environs, network environ analysis, environment, indirect effects, input-output analysis, 
connectivity, food web, trophic dynamics 



1. Introduction 

The dominance of indirect effects hypothesis is central to our new understanding of ecosystem ecology re- 
vealed through systems analyses like Ecological Network Analysis (ENA) (Ulanowicz, 1986; Fath and Patten, 
1999b; Dunne et al., 2002; Belgrano et al., 2005; Bascompte and Jordano, 2007; J0rgensen et al., 2007; McRae et al. 
2008; Christian et al., 2009; Baird et al., 2011). While empirical results arc also highlighting the important 
roles of indirect effects (e.g. Wootton, 1994; Berger et al., 2008; Dickottcr et al., 2007; Menendez et al., 
2007; Letnic et al., 2009; Walsh and Reznick, 2010), network analysis shows how the integral of indirect 
interactions across a whole web of interactions can transform the effective relationships between species 
so that they tend to be more positive (Ulanowicz and Puccia, 1990; Patten, 1991; Fath and Patten, 1998; 
Bondavalli and Ulanowicz, 1999), and more evenly distribute the system resources (Fath and Patten, 1999a; 
Borrett and Salas, 2010) and control (Schramski et al, 2006, 2007). 

When applying ENA, ecologists use the ratio of indirect to direct flows (I/D) to indicate the importance 
of indirect effects (Higashi and Patten, 1989; Fath and Patten, 1999b; J0rgensen et al., 2007). For example, 
Borrett et al. (2006) used I/D to show the temporal consistency of the organization of the Neuse River 
Estuary ecosystem, and Baird et al. (2009) found that indirect effects were dominant in the Sylt-R0mo Bight 
ecosystem regardless of the model aggregation scheme. Borrett and Freeze (2011) distinguished between two 
ways of calculating this ratio that they termed the unit and realized formulations. The unit calculation is the 
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traditional approach that is performed after a classic systems analysis technique is applied. In this case, each 
species is assumed to receive a single unit of boundary flux, and then the direct and indirect flow intensities 
in the system are summed. This unit calculation shows what the relative indirect flows would be if every 
node received the same input. It is useful because it facilitates within system comparisons (Whipple et al., 
2007) and due to its close relationship to the eigenvalues and vectors of underlying matrices (Borrett et al., 
2010). Alternatively, the realized calculation uses the observed boundary flow vector instead of an input 
vector of ones. This approach was first introduced by Borrett et al. (2006) but was clarified and formalized 
in Borrett and Freeze (2011). Its analytical advantage is that the system flows are scaled with the observed, 
typically non-homogeneous boundary fluxes. 

Borrett and Freeze (2011) showed that the quantitative and qualitative results could be changed by this 
shift in perspective, highlighting the importance of the system's connection to its broader environment. For 
example, the wet season model of carbon flux in the Everglades graminoid marshes (Ulanowicz et al., 2000; 
Hcymans et al., 2002) exhibits the dominance of indirect effects when analyzed with the realized method 
(I/D = 1.2) but does not when analyzed with the unit approach (I/D= 0.81). In this example, both the 
numerical value of the indicator and its interpretation changed. 

An important aspect of ENA throughflow analysis is that it has two orientations: input and output. 
Conceptually, we can either pull the energy-matter out of the system and trace its origin back through the 
system to the inputs, or we can push the inputs into the system and then trace where they travel through 
the system until they exit the system. The first type of analysis is termed the input analysis because in 
Leontief's (1966) original application he was determining what raw materials were required as inputs into 
the economic system to generate a single output, like a car. The second analysis is forward looking in 
direction, but it is termed the output analysis because we are following the inputs to their output from 
the system. In this case, it is the output generated that is the focus. Patten exploited this bidirectional 
analytical feature in the development of his environ concept (Patten, 1978, 1981, 1982) and the subsequent 
environmental theory. While some ENA components leverage the differences between the input and output 
orientation (e.g. Patten and Auble, 1981; Ulanowicz and Puccia, 1990; Fath and Patten, 1998; Patten, 1991; 
Schramski et al., 2006, 2007), the alternative perspectives can generate conflicting and confusing results. 

In this short communication we show that while the input and output oriented unit I/D can differ in 
quantitative and qualitative nature, the realized input and output oriented I/D are always identical. We 
first show this result numerically for a set of network ecosystem models drawn from the literature. We then 
provide a mathematical proof for why this relationship will generally hold when the models are at steady 
state. We conclude this paper by arguing that this identity bolsters the utility of the realized metric for (1) 
testing the generality of network hypotheses like the dominance of indirect effects (Salas and Borrett, 2011), 
(2) applying the ecological insights of this type of ENA to the original system, and (3) comparing ecosystems' 
organization (Baird et al., 1991, 2009; Borrett et al., 2006; Borrett and Osidele, 2007; Whipple et al, 2007). 

2. Ecological Network Analysis 

ENA is well described in the literature (e.g., Patten et al., 1976; Ulanowicz, 1986, 2004; Fath and Patten, 
1999b; Fath and Borrett, 2006; Schramski et al., 2010), but we recount the input and output throughflow 
analyses here as this is essential for our discussion. Borrett and Freeze (2011) focused on the distinction 
between the idealized and realized forms of the I/D ratio. Here we focus on contrasting the input and output 
oriented analyses. 

2.1. Model Input 

Ecologists apply ENA to network models of energy-matter storage and flux in ecosystems. These models 
are alternatively referred to as compartment models or energy-matter budgets. In the general form of these 
models, n nodes represent species, groups of species, or abiotic components, and the L weighted directed 
edges represent the flow of energy-matter generated by some ecological process (e.g., consumption, excretion, 
harvesting). Let F nx „ = (fij) represent the observed flow from ecosystem compartment j to compartment 
i (e.g., j —> i), z n xi be a vector of node inputs originating from outside the system, and yi xn be a vector 
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of flows from each node that exit the system. Sometimes the outputs are subdivided into those that are 
thermo dynamically degraded and those that reflect high quality energy-matter exported from the system 
(e.g., Ulanowicz, 1986). For the analysis presented here, we lump these losses together because once they 
cross the system boundary our analysis is blind to their fate, making this distinction less critical. 

ENA throughflow analysis assumes that the model traces a single conserved currency (e.g., energy, 
nitrogen, phosphorus), and that the system is at a steady state (inputs equal outputs). The steady state 
assumption is crucial for the throughflow path decomposition (Borrett et al., 2010). 

For this research, we examined 50 ecosystem network models that represent 35 distinct primarily fresh- 
water and marine ecosystems (Table 1). The models exhibit a range of sizes (4 < n < 125), connectance 
(number of direct links divided by the total possible number of links; 0.03 < (C = L/n 2 ) < 0.40), and 
recycling (Finn Cycling Index (Finn, 1980), which is the proportion of total system throughflow derived 
from recycled flux; < FCI < 0.51). The core of these trophically-based networks is a food web, but 
other ecological processes are included and model compartments may represent non-living resource pools 
such as particulate organic carbon. We applied Allesina and Bondavalli's (2003) AVG2 balancing algo- 
rithm to 15 models that were not initially at steady state because this technique tends to generate the 
least distortion of network properties. This is the same set of models used in several recent network stud- 
ies (Borrett and Salas, 2010; Borrett and Freeze, 2011; Salas and Borrett, 2011), and the model data are 
available from http://people.uncw.edu/borretts/research.html. 

2.2. Throughflow Analysis 

Input and output oriented throughflow analysis have three main steps. First, we determine the through- 
flow vector T, which is the total amount of energy-matter flowing into or out of each node. This can be 
calculated from the initial model information as follows: 

n 

7f = E^+ z > (i = l,2,...,n),and (1) 
i=i 

n 

IT £/y ■ (J = l,2,...,n). (2) 

i=l 

At steady state, T { U1 = (T^ out ) T = T nx i = (Tj). From this vector, we derive the first whole-system indicator, 
total system throughflow (TST = J2j=i Tj). TST indicates the total magnitude of flow activity. 

Second, we calculate the input G' rlX n = {g[j) and output G nxn = (cjij) direct flow intensities from node 
j to i. These are defined as 

gfa ee / y /2f , and (3) 
g tJ ee fa/T?^. (4) 

Here, g'^ is the fraction of input at receiver node i contributed from the donor node j, while gij is the 
fraction of output throughflow at donor node j contributed to node i. The subtle distinction between the 
input and output oriented direct flow intensities in equations (3 and 4) is whether the observed flow is 
normalized by the throughflow of the receiving node (input case) or the donating node (output case). The 
g'ij and gij values are dimensionless and the column sums of the matrices must lie between and 1 because 
of thermodynamic constraints of the original model (see Borrett et al., 2010, for details). 
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Third, we determine the integral flow intensities N' = (nij) (input) and N = ijiij) (output) as 

■oo 

N'=^G""= ^ + G^+ G ,2 + ... + G' m + ... , (5) 

Boundary Direct Indirect 

and 

oo 

N=^G m = I + G^+ G 2 + ... + G m + ... . (6) 

m =0 Boundary Direct Indirect 

In equations (5 and 6) , I = (iy ) is the matrix multiplicative identity and the elements of G' m and G m are 
the fractions of boundary flow that travels from node j to i over all pathways of length m. The exact values 
of N' and N can be found using the following identities because given our model definitions the power series 
must converge. Thus, 

N' = (I- G') -1 , and (7) 
N = (I-G) _1 . (8) 

The n^j and n,j elements represent the intensity of boundary flow that passes from j to i over all pathways 
of all lengths. These values integrate the boundary, direct, and indirect flows. 

The differences between N' and N results from the input and output perspectives. However, we can use 
either of these integral flow matrices to recover T as follows: 

f = yN', (9) 
f = Nz. (10) 

2.3. Unit and Realized Indirect Effects 

There are two methods of calculating I/D used in ENA to quantify the importance of indirect flows. We 
refer to these as the unit and realized methods. The unit method assumes that each node receives a single 
unit of input, which we will represent as a vector of ones with length equal to the number of nodes [l] nx i 
or output [l]ixn- The unit input and output I/D metrics are calculated as follows. 

-'/-k'unit, input — ~ T= : l-LJ-J 

EILi([i]ix„G') 4 

Tin ^((N-i-GitLi), ri2 . 

l/Uxmlt, output - — " [U) 

In contrast, the realized method uses the observed boundary vectors z or y. 

I/D SUflN'-I-G')), 

1 1 ^realized, input — rnn ( -r~i i\ K 10 ) 

£" ((N - 1 - G)z) . 

1 / ^realized, output — -#\ V / 

As Borrett and Freeze (2011) discuss, the magnitude of the boundary flow vectors, be it a unit vector 
[l] nxl or the observed boundary vectors z or y, cancels in the ratio measure. However, the different 
distribution of boundary flows differentially excites the flow intensities in the integral flow matricics. It is 
the distribution of the boundary inputs and outputs that is key. 
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3. Results 



We present our results in two parts. The first section shows the numerical results of applying EN A and 
the alternative formulations of I/D to the 50 network models. The second section summarizes the analytical 
work that generalizes this result. 

3.1. Numerical Results 

Figure 1 compares the input and output oriented I/D values of the network models. Panel (a) illustrates 
the variation that can occur between the input and output oriented I/D values when using the idealized 
formulation. The adjusted r 2 value for a linear regression fit to this data set is 0.82, and the Pearson's 
correlation coefficient is 0.91. To ensure our balancing routine did not heavily influence our results, we also 
calculated the Pearson's correlation coefficient for only the 35 models not balanced and found it to be 0.92. 
Panel (b) shows the one to one correspondence between the input and output I/D values when we use the 
realized formulation. These values are identical and have a Pearson's correlation coefficient of 1. 

3.2. Analytical Results 

The steady state assumption actually implies very strong compatiblity conditions on the powers of G 
and G', namely that for each nonnegative integer m, the sizes of G m z and yG m must coincide in the sense 
that we have the equality 

n n 

£(G" l z%- = £(yG'™) i , (15) 

i=l i=l 

as demonstrated in the proof of Lemma 6.3 in Appendix A. 

The terms of the numerator and denominator of realized I/D each are of the form found on the left hand 
side in equation (15) for realized output I/D or of the form found on the right hand side in equation (15) 
for realized input I/D . It follows that at steady state, the realized input I/D and output I/D must be 
identical. This mathematical generalization of the numerically observed results is recorded in Theorem 6.4 
of Appendix A. 

4. Discussion 

In this note we show that using the realized formulation the input and output oriented I/D values 
are identical when the system is at steady state. Like Bata et al. (2007), this work provides a conceptual 
and mathematical condensation for Ecological Network Analysis. This is part of a methodological matu- 
ration process in which investigators are finding relationships amongst the alternative analyses (input vs. 
output or throughflow vs. storage) and reconsidering initial assumptions that have become canalized (e.g. 
Schramski et al., 2010; Matamba et al., 2009). Here, we are not fully pruning the analytical options, but 
clarifying the consequences of alternative methods. When an investigator selects to use the realized method 
as described by Borrett and Freeze (2011), one advantage is that the input and output values are identical. 

The realized formulation of I/D functions as one description of the internal organization of the ecosys- 
tem. It is intuitively satisfying that our understanding of this organization does not change because our 
point of view switches (input or output). Further, as the realized formulation considers both the internal 
system organization and the external environmental connection, it characterizes the complete system as 
it was initially described by the model builders. Thus, the realized metric is most appropriate for test- 
ing the generality of systems ecology hypotheses like the dominance of indirect effects (Salas and Borrett, 
2011), applying ENA insights back to the original system (Zhang et al., 2010a,b), and comparing ecosystem 
organization (Baird et al., 1991; Whipple et al., 2007; Ray, 2008). 

A limitation of this work is that the realized input and output I/D should only be identical when the 
observed system is at steady state, as shown in the mathematical proof. This means systems not at steady 
state will likely have different relative contributions of indirect flows when observed forward through the 
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network instead of backwards through the network. This is evident when we recall equations (9) and (10) 
and that T m = T out at steady state. Without this identity (i.e., not at steady state), there is no reason 
to expect the mathematical decomposition of the input and output throughflow into boundary, direct, and 
indirect flows to be identical. We suspect that network particle tracking as implemented in EcoNet (Kazanci, 
2007; Tollner et al., 2009) will be useful for further exploring the consequences of these input and output 
perspective differences. 

This work is one more essential step in the maturation of the ecosystem and environmental theory devel- 
oping in association with Ecosystem Network Analysis (e.g. Abarca- Arenas and Ulanowicz, 2002; Allesina and Ulanowicz, 
2004; Borrctt et al., 2010; Kaufman and Borrett, 2010; Schramski et al., 2010). Multiple lines of investiga- 
tion are building confidence in ENA results including the initial empirical validation by Dame and Christian 
(2008) and the agreement between the Eulerian approach used here and the Lagrangian approach of net- 
work particle tracking (Matamba et al., 2009), but continued rigorous empirical testing of network ecology 
hypotheses is essential. Loehle (1987) points out that this theory maturation is essential so that theory 
predictions clarify to the point that robust empirical test be successful. 
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6. Appendix A 

The steady state assumption that f in = (f out ) T means that Tf 1 = T £ out for all indices £. We have the 
following proposition as an immediate consequence. 

Proposition 6.1. When the observed system is at steady state, we have that y]AGz)j = ^2i(yG')i. 
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Proof. Note that 



E( G ^ = E 

i i 

= E 



E 



J nnout 



J mi] 



e[(V E(/»)) 
E(/«) - E f 



(by equation (1)) 



fijfjk 



E(/«)-E 



Efo) 



= E(/«)-E 



3, ft 



fjk 



E(/«)-E(/i*)+E( 

*J J.ft 3>ft V 



(by equation (2)) 



^ ; I <TXmt % 

j,ft V i 



E(^ G ')* 



We extend the observation of Proposition 6.1. 
Lemma 6.2. For each integer m > 2, we have that 



and 



(yG™),= £ 



11 I rpOUt 



"pout 
e i n=l 



£ tn—X/p 



„=i \ 



Proof. We observed in the proof of Proposition 6.1 that 



J rpOUt 



□ 
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It follows that 



(G 2 z) 4 = (G(G2))< 



E 
E 
E 



E 



fij fjk \ 

^OUt J^QUt j 
j k J 

fit! flll2 



and by iteration we obtain the stated formula for (G m z),. The formula for (yG m \ is similarly obtained. □ 

The following lemma makes use of Lemma 6.2 to extend the observations of Proposition 6.1. 
Lemma 6.3. For each integer m > 0, we have that 



i=i 



Proof. The case m = 0, namely that = j/j, follows directly from the steady-state assumption, and 

the case m = 1 was observed in Proposition 6.1. Suppose then that m > 2. 

Observe that we have the following identity by re-indexing the left-hand side via it >-> tt+i for 1 < t < 
in — 1, £ m t-> £ x , followed by renaming £x to l m +\ and then i to i\\ 



E 



m — 1 



= E 



(16) 
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From Lemma 6.2 it follows that 



(recall that T m = T out ) 




n 



/<„<„+! 



(relabelling i m +i as i) 



□ 



Theorem 6.4. Realized input and output I/D are identical when the observed system is at steady state. 



Proof. Lemma 6.3 implies the equality of the denominators of realized input and output I/D as well as 
equality of corresponding summands of the numerators of realized input and output I/D . □ 
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Figure 1: Relationship between input and output oriented measures of the ratio of indirect-to-direct effects calculated using 
the (a) unit and (b) realized metrics. 
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Table 1: Fifty trophically-bascd ecosystem network models. 



Model 


units 




nt 


ct 


TSTt 




Source 


Lake Findley 


gC vcT 2 


yr- 1 


■1 


0.38 


51 


0.30 


Richey et al. (1978) 


Mirror Lake 


gC m -2 


yr- 1 


5 


0.36 


218 


0.32 


Richey et al. (1978) 


Lake Wingra 


gC m~ 2 


yr- 1 


5 


0.40 


1,517 


0.40 


Richey et al. (1978) 


Marion Lake 


gC m™ 2 


yr" 1 


5 


0.36 


243 


0.31 


Richey et al. (1978) 


Cone Springs 


kcal m - 


2 yr" 1 


5 


0.32 


30,626 


0.09 


Tilly (1968) 


Silver Springs 


kcal m - 


2 yr^ 1 


5 


0.28 


29,175 


0.00 


Odum (1957) 


English Channel 


kcal m — 


2 yi^ 1 


6 


0.25 


2,280 


0.00 


Brylinsky (1972) 


Oyster Reef 


kcal m - 


2 yr -i 


6 


0.33 


84 


0.11 


Dame and Patten (1981) 


Sommc Estuary 


mgC m" 


-2 d -i 


9 


0.30 


2,035 


0.14 


Rybarczyk et al. (2003) 


Bothnian Bay 


gC m" 2 


yr" 1 


12 


0.22 


130 


0.18 


Sandberg et al. (2000) 


Bothnian Sea 


gC m -2 


yr" 1 


12 


0.24 


458 


0.27 


Sandberg et al. (2000) 


Ythan Estuary 


gC m~ 2 


yr- 1 


13 


0.23 


4,181 


0.24 


Baird and Milne (1981) 


Baltic Sea 


mgC rrT 


' 2 d _1 


15 


0.17 


1,974 


0.13 


Baird et al. (1991) 


Ems Estuary 


mgC m" 


-2 d -l 


15 


0.19 


1,019 


0.32 


Baird et al. (1991) 


Swarkops Estuary 


mgC m" 


-2 d -l 


15 


0.17 


13,996 


0.47 


Baird ct al. (1991) 


Southern Benguela Upwclling 


mgC m" 


-2 d -l 


16 


0.23 


1,774 


0.19 


Baird ct al. (1991) 


Peruvian Upwelling 


mgC rrT 


-2 d -l 


16 


0.22 


33,496 


0.04 


Baird ct al. (1991) 


Crystal River (control) 


mgC m~ 


-2 d -l 


21 


0.19 


15,063 


0.07 


Ulanowicz (1986) 


Crystal River (thermal) 


mgC m~ 


-2 d -l 


21 


0.14 


12,032 


0.09 


Ulanowicz (1986) 


Charca de Maspalomas Lagoon 


mgC m~ 


-2 d -l 


21 


0.13 


6,010,331 


0.18 


Almunia et al. (1999) 


Northern Benguela Upwclling 


mgC m" 


-2 d -l 


24 


0.21 


6,608 


0.05 


Heymans and Baird (2000) 


Neuse Estuary (early summer 1997) 


mgC m~ 


' 2 d" 1 


30 


0.09 


13,826 


0.12 


Baird et al. (2004b) 


Neuse Estuary (late summer 1997) 


mgC m~ 


' 2 d" 1 


30 


0.11 


13,038 


0.13 


Baird ct al. (2004b) 


Neuse Estuary (early summer 1998) 


mgC m~ 


-2 d -i 


30 


0.09 


14,025 


0.12 


Baird et al. (2004b) 


Neuse Estuary (late summer 1998) 


mgC m~ 


-2 d ~l 


30 


0.10 


15,031 


0.11 


Baird et al. (2004b) 


Gulf of Maine 


g ww m 


-2 yr -i 


31 


0.35 


18,382 


0.15 


Link ct al. (2008) 


Georges Bank 


g ww m 


- 2 yr- 1 


31 


0.35 


16,890 


0.18 


Link ct al. (2008) 


Middle Atlantic Bight 


g ww m 




32 


0.37 


17,917 


0.18 


Link et al. (2008) 


Narragansett Bay 


mgC m" 


-2 yr" 1 


32 


0.15 


3,917,246 


0.51 


Monaco and Ulanowicz (1997) 


Southern New England Bight 


g ww m 


-2 yr" 1 


33, 


0.03 


17,597 


0.16 


Link et al. (2008) 


Chesapeake Bay 


mgC m" 


-2 yr" 1 


36 


0.09 


3,227,453 


0.19 


Baird and Ulanowicz (1989) 


St. Marks Seagrass, site 1 (Jan) 


mgC rrT 


-2 d -l 


51 


0.08 


1,316 


0.13 


Baird ct al. (1998) 


St. Marks Seagrass, site 1 (leb) 


mgC rrT 


-2 J —1 


51 


0.08 


1,591 


0.11 


Baird et al. (1998) 


St. Marks Seagrass, site 2 (Jan) 


mgC m" 


-2 d -l 


51 


0.07 


1,383 


0.09 


Baird et al. (1998) 


St. Marks Seagrass, site 2 (Feb) 


mgC m~ 


-2 d -l 


51 


0.08 


1,921 


0.08 


Baird ct al. (1998) 


St. Marks Seagrass, site 3 (Jan) 


mgC m" 


-2 d -l 


51 


0.05 


12,651 


0.01 


Baird ct al. (1998) 


St. Marks Seagrass, site 4 (Feb) 


mgC rrT 


-2 d -l 


51 


0.08 


2,865 


0.04 


Baird et al. (1998) 


Sylt R0m0 Bight 


mgC m~ 


-2 d -l 


59 


0.08 


1,353,406 


0.09 


Baird et al. (2004a) 


Graminoids (wet) 


gC m" 2 


yr -1 


66 


0.18 


13,677 


0.02 


Ulanowicz et al. (2000) 


Graminoids (dry) 


gC m- 2 


yr -1 


66 


0.18 


7,520 


0.04 


Ulanowicz et al. (2000) 


Cypress (wet) 


gC m~ 2 


yr" 1 


68 


0.12 


2,572 


0.04 


Ulanowicz et al. (1997) 


Cypress (dry) 


gC m- 2 


yr- 1 


68 


0.12 


1,918 


0.04 


Ulanowicz et al. (1997) 


Lake Oneida (prc-ZM) 


gC m -2 


yr- 1 


71 


0.22 


1,638 


< 0.01 


Miehls et al. (2009a) 


Lake Quintc (pre-ZM) 


gC m -2 


yr" 1 


71 


0.21 


1,467 


< 0.01 


Miehls et al. (2009b) 


Lake Oneida (post-ZM) 


gC m" 2 


yr" 1 


76 


0.22 


1,365 


< 0.01 


Miehls et al. (2009a) 


Lake Quinte (post-ZM) 


gC m -2 


yr" 1 


80 


0.21 


1,925 


0.01 


Miehls et al. (2009b) 


Mangroves (wet) 


gC m- 2 


yr- 1 


94 


0.15 


3,272 


0.10 


Ulanowicz et al. (1999) 


Mangroves (dry) 


gC m -2 


yr— 1 


91 


0.15 


3,266 


0.10 


Ulanowicz et al. (1999) 


Florida Bay (wet) 


mgC m~ 


" 2 yr" 1 


125 


0.12 


2,721 


0.14 


Ulanowicz et al. (1998) 


Florida Bay (dry) 


mgC m" 


- 2 yr" 1 


125 


0.13 


1,779 


0.08 


Ulanowicz et al. (1998) 



' n is the number of nodes in the network mode 1, C = L/n 2 is the model conncctance when L is the number of direct links or 
energy— matter transfers, TST — S fij z i ls ^ ne t°t a l system throughflow, and FCI is the Finn Cycling Index (Finn, 1980). 
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